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The sonic analog of a gravitational black hole in dilute-gas Bose-Einstein condensates is investigated. 
It is shown that there exist both dynamically stable and unstable configurations which, in the 
hydrodynamic limit, exhibit a behavior completely analogous to that of gravitational black holes. 
The dynamical instabilities involve creation of quasiparticle pairs in positive and negative energy 
states. We illustrate these features in two qualitatively different one-dimensional models, namely, 
a long, thin condensate with an outcoupler laser beam providing an "atom sink," and a tight ring- 
shaped condensate. We have also simulated the creation of a stable sonic black hole by solving 
the Gross-Pitaevskii equation numerically for a condensate subject to a trapping potential which 
is adiabatically deformed. A sonic black hole could in this way be created experimentally with 
state-of-the-art or planned technology. 
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I. INTRODUCTION 

Many investigations of dilute gas Bose-Einstein con- 
densates arc directed towards experimentally creating 
nontrivial configurations of the semiclassical mean field, 
or to predicting the properties of such configurations in 
the presence of quantum fluctuations. Such problems are 
hardly peculiar to condensates: the quantum neighbor- 
hoods of interesting classical backgrounds are important 
areas of research in most fields of physics. But ultracold 
dilute gases are so easy to manipulate and control, both 
experimentally Q| and theoretically ||], that they may 
allow us to decipher less amenable systems by analogy. 
As an essay in such an application of condensates, in this 
paper we discuss the theoretical framework and propose 
an experiment to create the analog of a black hole in the 
laboratory and simulate its radiative instabilities. 

It is now commonly believed that, even in the con- 
text of elementary particle physics, quantum field the- 
ory arises from a still unknown underlying structure: it 
is an effective dynamical theory, describing the low en- 
ergy limit of collective phenomena of the underlying mi- 
croscopic theory. From this viewpoint, our description 
of (more) fundamental phenomena, such as gravity or 
electromagnetism, is actually similar to the theoretical 
descriptions of many phenomena of condensed matter. 
To understand superfluidity, superconductivity or dilute 
Bose-Einstein condensation, we describe the dynamics 
of the system in terms of collective modes (quasipar- 
ticles) whose typical size is much larger than the dis- 
tances between the particles that constitute the under- 
lying medium; but even electrons and photons must be 
considered as the quasiparticles of the deeper theory we 
do not yet know. In this sense we may say that the 
major difference between our fundamental theories and 



those we use in condensed matter is that in the latter 
case the next microscopic level of description is actually 
well understood. 

With this fundamental background in mind, it is not 
so surprising that condensed matter analogs of nontriv- 
ial configurations appearing in relativistic quantum field 
theories and gravitation can be constructed. For ex- 
ample, 3 He has been proposed as a laboratory counter- 
part of high-energy particle physics. It has been argued 
that, under appropriate conditions, excitations around 
the ground state of the system may resemble the particle- 
spectrum of gauge theories of high energy physics ||. 
These condensed matter systems have also been used to 
simulate topological defects characteristic of gauge theo- 
ries and which are considered to have played a cosmologi- 
cal role in the early stages of the evolution of the universe 
such as monopoles and cosmic strings ||. 

The past decade has witnessed an increasing interest 
in simulating gravitational configurations and processes 
in condensed matter systems in the laboratory. The key 
observation was originally made by Unruh and fur- 
ther analyzed by Visser j(||7j : phononic propagation in a 
fluid is described by a wave equation which, under ap- 
propriate conditions, can be interpreted as propagation 
in an effective relativistic curved spacetime background, 
the spacetime metric being entirely determined by the 
physical properties of the fluid under study, namely, its 
density and flow velocity. Unruh urged a specific moti- 
vation |Q| for examining the hydrodynamic analogue of 
an event horizon namely that as an experimentally 
and theoretically accessible phenomenon it might shed 
some light on the Hawking effect |j) (thermal radiation 
from black holes, stationary insofar as the back reaction 
is negligible). In particular, one would like to gain in- 
sight into the role in the Hawking process of ultrahigh 



1 



Garay, Anglin, Cirac, and Zoller. Sonic black holes in dilute Bose-Einstein condensates 



2 



frequencies |,0[T|]. 

An event horizon for sound waves appears in princi- 
ple wherever there is a surface through which a fluid 
flows at the speed of sound, the flow being subsonic 
on one side of the surface and supersonic on the other. 
There is a close analogy between sound propagation on 
a background hydrodynamic flow, and field propagation 
in a curved spacetime; and although hydrodynamics is 
only a long-wavelength effective theory for physical (su- 
per)ffuids, so also field theory in curved spacetime is to 
be considered a long-wavelength approximation to quan- 
tum gravity |5|,|^]- Determining whether and how sonic 
black holes radiate sound, in a full calculation beyond the 
hydrodynamic approximation or in an actual experiment, 
can thus offer some suggestions about black hole radiance 
and its sensitivity to high frequency physics (beyond the 
Planck scale). The possibility that such high frequencies 
might have consequences for observably low frequency 
phenomena is one of the main reasons that black holes 
have deserved much attention: there is reason to expect 
that an event horizon can act as a microscope, giving us a 
view into physics on scales below the Planck length. This 
is because modes coming from an event horizon are red- 
shifted into the low-energy regime as they propagate out 
to be observed far away from the black hole. Conversely, 
if we imagine tracking the observed signal back towards 
its source, the closer we get to the horizon, the shorter 
the wave length of the signal must be, until at the very 
horizon we must either reach infinite energy scales, or en- 
counter a breakdown in general relativity and quantum 



field theory in curved spacetime |12|. 

In understanding this problem, hydrodynamic and 
condensed matter analogs of black holes may offer some 
of the experimental guidance otherwise difficult to obtain 
in the case of gravity fl^l . Under appropriate conditions 
and approximations (which can be basically summarized 
in the requirement the wave lengths of the perturbations 
be sufficiently large), the propagation of collective fluc- 
tuations (phonons) admits an effective general relativis- 
tic description, in terms of a spacetime metric. This 
long wavelength regime would correspond analogically to 
quantum field theory in curved spacetime. The effective 
phonon metric may describe black holes, as in general rel- 
ativity, and so a phonon Hawking effect may be possible; 
and certainly the problem of arbitrarily high frequencies 
at the horizon is also present. But in this case, when at 
short wave lengths the metric approximation is no longer 
valid and a more microscopic theory must be used in- 
stead, the accurate microscopic theory is actually known. 
And if the hydrodynamic system is a dilute Bose-Einstein 
condensate, the microscopic theory is actually tractable 
enough that we can make reliable calculations from first 
principles. As we will argue, trapped bosons at ultralow 
temperature can indeed provide an analogue to a black- 
hole spacetime. Similar analogues have been proposed in 
other contexts, such as superfluid helium [ff3|, solid state 



physics 1 14 1, and optics |ll|; but the outstanding recent 
experimental progress in cooling, manipulating and con- 
trolling atoms |l6| make Bose-Einstein condensates an 
especially powerful tool for this kind of investigation. 

The basic challenge of our proposal is to keep the 
trapped Bose-Einstein gas sufficiently cold and well iso- 
lated to maintain a locally supersonic flow long enough 
to observe its intrinsic dynamics. Detecting thermal 
phonons radiating from the horizons would obviously be 
a difficult additional problem, since such radiation would 
be indistinguishable from many other possible heating 
effects. This further difficulty does not arise in our pro- 
posal, however, because the black-hole radiation we pre- 
dict is not quasistationary, but grows exponentially under 
appropriate conditions. It should therefore be observable 
in the next generation of atom traps, and may also raise 
new issues in the theory of gravitational black holes. 

In this paper, we extend and generalize the results of 
Rcf . , including a more detailed analysis of the model 
considered there as well as a new qualitatively different 
case. The paper is organized as follows. In Sec. II, we 
show how sonic horizons in dilute Bose-Einstein conden- 
sates may appear in the hydrodynamic approximation, 
discuss the regime of validity of such approximation, as 
well the validity of one-dimensional models. Section III is 
devoted to the study of sonic horizons in condensates sub- 
ject to tight ring-shaped external potentials. We present 
numerical results showing that both stable and unstable 
black holes may be created under realistically attainable 
conditions in current or near- future laboratories. We also 
study the nature of the dynamical instabilities that ap- 
pear for certain configurations. In Sec. IV, we discuss a 
different configuration, namely that of a sink-generated 
black hole in an infinite one-dimensional condensate, and 
show that there also exist black-hole configurations, al- 
though they are not stable. We summarize and conclude 
in Sec. V. The Appendix is devoted to the issues of 
redundancy and normalization of the dynamically unsta- 
ble Bogoliubov modes (associated with complex eigenfre- 
quencies) . 

II. SONIC BLACK HOLES IN CONDENSATES 

A Bose-Einstein condensate is the ground state of a 
second quantized many body Hamiltonian for N inter- 
acting bosons trapped by an external potential V4 x t(x) 
At zero temperature, when the number of atoms is 
large and the atomic interactions are sufficiently small, 
almost all the atoms are in the same single-particle quan- 
tum state ^(x, t), even if the system is slightly perturbed. 
The evolution of '5 is then given by the well-known Gross- 
Pitaevskii equation, which in appropriate units can be 
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written as 

\ 2m to ) 

where to is the mass of the individual atoms and a is 
the scattering length. The wave function of the con- 
densate is normalized to the total number of atoms 
/d 3 x|*(x,i)| 2 = TV. 

Our purposes do not require solving the Gross- 
Pitaevskii equation with some given external potential 
V ex t (x) ; our concern is the propagation of small collective 
perturbations of the condensate, around a background 
stationary state 

* s (x,t) = ^wWe-W, 

where \x is the chemical potential. Thus it is only nec- 
essary that it be possible, in any external potential that 
can be generated, to create a condensate in this state. 
Indeed, many realistic techniques for "quantum state en- 
gineering," to create designer potentials and bring con- 
densates into specific states, have been proposed, and 
even implemented successfully |l6| ; and our simulations 
indicate that currently known techniques should suffice 
to generate the condensate states that we propose. 

Perturbations about the stationary state ^> s (x, t) obey 
the Bogoliubov system of two coupled second order dif- 
ferential equations. Within the regime of validity of the 
hydrodynamic (Thomas- Fermi) approximation J2|, these 
two equations for the density perturbation g and the 
phase perturbation <p in terms of the local speed of sound 



c(x) = — ^47rap(x), 

TO 

and the background stationary velocity 
v= -Vtf 

TO 

read 

_ / to 9 \ • Airah 

g=-V- - — -C 2 V(j) + VQ) = -v-V0 g. 

\AiTan / to 

Furthermore, low frequency perturbations are essentially 
just waves of (zero) sound. Indeed, the Bogoliubov equa- 
tions may be reduced to a single second order equation 
for the condensate phase perturbation <fi. This differen- 
tial equation has the form of a relativistic wave equation 
dfii^f—gg^ducj)) = 0, with g = detg M „, in an effective 
curved spacetime with the metric g^ v being entirely de- 
termined by the local speed of sound c and the back- 
ground stationary velocity v. Up to a conformal factor, 
this effective metric has the form 




This class of metrics can possess event horizons. For 
instance, if an effective sink for atoms is generated at the 
center of a spherical trap (such as by an atom laser out- 
coupling technique |l^] ) , and if the radial potential profile 
is suitably arranged, we can produce densities p(r) and 
flow velocities v(x) = —v(r)r/r such that the quantity 
c 2 — v 2 vanishes at a radius r — Th, being negative inside 
and positive outside. The sphere at radius is a sonic 
event horizon completely analogous to those appearing 
in general relativistic black holes, in the sense that sonic 
perturbations cannot propagate through this surface in 
the outward direction [0J|,[jJ . This can be seen explicitly 
by writing the equation for the radial null geodesies of 
the metric g^ v : 

r± = —v ± c, 

which can be obtained from setting the proper interval 
ds 2 = g^dx^dx" equal to zero and restricting the al- 
lowed motion to the radial direction, so that 

-(c 2 -v 2 ) + 2vf + f 2 =0. 

The ingoing null geodesic r_(t) is not affected by the 
presence of the horizon and can cross it in a finite coor- 
dinate time t. The outgoing null geodesic r + (t) on the 
other hand needs an infinite amount of time to leave the 
horizon since r+ = at the horizon. The physical mech- 
anism of the sonic black hole is quite simple: inside the 
horizon, the background flow speed v is larger than the 
local speed of sound c, and so sound waves are inexorably 
dragged inwards. 

In fact there are two conditions which must hold for 
this dragged sound picture to be accurate. Wavelengths 
larger than the black hole itself will of course not be 
dragged in, but merely diffracted around it. And pertur- 
bations must have wavelengths 

A> — , ; 

mc mc-y/\l-v/c\ 

Otherwise they do not behave as sound waves since they 
lie outside the regime of validity of the hydrodynamic 
approximation. These short-wavelength modes must be 
described by the full Bogoliubov equations, which allow 
signals to propagate faster than the local sound speed, 
and thus permit escape from sonic black holes. So, to 
identify a condensate state \P S as a sonic black hole, there 
must exist modes with wavelengths larger than these 
lower limits (which in terms of the local healing length 
£(x) = ?L/[toc(x)] read A > 2tt£, 2tt£/- v /|1 - v/c\), but 
also smaller than the black hole size. Even if such an 
intermediate range does exist, the modes outside it may 
still affect the stability of the black hole as discussed be- 
low. 

As it stands, this description is incomplete. The con- 
densate flows continually inwards and therefore at r = 
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there must be a sink that takes atoms out of the con- 
densate. Otherwise, the continuity equation V(/ov) = 0, 
which must hold for stationary configurations will be vi- 
olated. From the physical point of view, such a sink can 
be accomplished by means of an outcoupler laser beam at 
the origin. (Such outcouplers are the basic mechanisms 
for making trapped condensates into "atom lasers," and 
they have already been demonstrated experimentally by 
several groups. A tightly focused laser pulse changes the 
internal state of the atoms at a particular point in the 
trap, and can also be made to give them a large momen- 
tum impulse. This ejects them so rapidly through the 
always dilute condensate cloud that they do not signifi- 
cantly disturb it; effectively, they simply disappear.) 

We have analyzed several specific systems which may 
be suitable theoretical models for future experiments, 
and have found that the qualitative behavior is analo- 
gous in all of them. Black holes which require atom sinks 
are both theoretically and experimentally more involved, 
however; moreover, maintaining a steady transonic flow 
into a sink may require either a very large condensate 
or some means of replenishment. We will therefore first 
discuss an alternative configuration which may be exper- 
imentally more accessible and whose description is par- 
ticularly simple: a condensate in a very thin ring that 
effectively behaves as a periodic one-dimensional system 
(Fig. |l|). Under conditions that we will discuss, the super- 
sonic region in a ring may be bounded by two horizons: 
a black hole horizon through which phonons cannot exit, 
and a 'white hole' horizon through which they cannot en- 
ter. Then we will analyze another simple one-dimensional 
model, of a long, straight condensate with an atom sink 
at the center (Fig. |^). 



White hole horizon 




Condensate cloud 

FIG. 1. The tight ring-shaped configuration, with both 
black and white horizons, and no singularity. Arrows indi- 
cate condensate flow velocity, with longer arrows for faster 
flow. 

The existence of instabilities that do not show up in the 
one-dimensional approximation is an important question 
in condensate physics, which is under active theoretical 
and experimental investigation. The essential principles 



have long been clear, inasmuch as the current dilute con- 
densates really are the weakly interacting Bose gases that 
have been used as toy models for superfluidity for several 
decades. The fact that actual critical velocities in liquid 
helium are generally far below the Landau critical veloc- 
ity is understood to be due partly to the roton feature of 
the helium dispersion relation, but this is not present in 
the dilute condensates. Viscosity also arises due to sur- 
face effects, however, and these may indeed afflict dilute 
condensates as well. The point here is that in addition to 
the bulk phonon modes considered by Landau, and quite 
adequately represented in our one-dimensional analysis, 
there may in principle be surface modes, with a different 
(and generally lower) dispersion curve. If such modes 
exist and are unstable, it is very often the case that, as 
they grow beyond the perturbative regime, they turn into 
quantized vortices, which can cut through the supercur- 
rent and so lower it. 




j 'Singularity' 



Biack hole horizons 



Outcoupled 
'atom laser' beam 



FIG. 2. The tight cigar-shaped configuration, with two 
black hole horizons and a 'singularity' where condensate is 
outcoupled. Arrows indicate condensate flow velocity, with 
longer arrows for faster flow. 

Whether or not such unstable surface modes actually 
exist in the Bogoliubov spectrum of a dilute condensate is 
an issue that has recently been analyzed both numerically 
and analytically, and it is quite clear that such surface 
modes exist only if the confining potential is quite rough 
(which is not only easy to avoid with a magnetic or opti- 
cal trapping field, but very hard to achieve) jl9) , or if the 
condensate dynamics in the directions perpendicular to 
the flow is hydrodynamic. That is, the condensate must 
be at least a few healing lengths thick, so that surface 
modes decaying on the healing length scale can satisfy 
all the required boundary conditions pQ] . By saying that 
we are considering an effectively one-dimensional conden- 
sate, we mean precisely that this is not the case. For 
instance, for the tight-ring model, in this regime, the ra- 
dial trap scale is the shortest length scale in the problem, 
and the radial trap frequency is the highest frequency; 
this effectively means that excitations of nontrivial radial 
modes, including surface modes, are energetically frozen 
out. (In the limit of radial confinement within the scat- 
tering length, our model breaks down for other reasons 
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- but the scattering length can easily be two orders of 
magnitude smaller than the healing length) . The issue of 
the supercurrent stability in tightly confined ring shaped 
traps has been addressed in Ref. |2l[] where the authors 
arrive at a positive conclusion and also clarify the role of 
finite temperature and possible trap anisotropy. 



on wavelengths shorter than the inner and the outer re- 
gions are indeed phononic implies U 3> 27r, which in turn 
requires w ^> 1 and 1 > 6 > 1/w 2 . In fact, detailed 
analysis shows that w ^ 5 is sufficient. 

A. Stability 



III. SONIC BLACK/WHITE HOLES IN A RING 

In a sufficiently tight ring-shaped external potential 
of radius R, motion in radial (r) and axial (z) cylin- 
drical coordinates is effectively frozen. We can then 
write the wave function as fy(z,r, 0,r) — f(z,r)&(0,r) 
and normalize $ to the number of atoms in the con- 
densate J Q T d8\<&(9)\ 2 = N, where with the azimuthal 
coordinate 9 we have introduced the dimensionless time 
t = (h/mR 2 )t. The Gross-Pitaevskii equation thus be- 
comes effectively one-dimensional: 



"VT>= [-\dl 



Ve 



"I 

N 



(i) 



where U = AnaNR 2 J dzdrr\f(z, r)\ 4 and V cxt (<9) is 
the dimensionless effective potential (in which we have 
already included the chemical potential) that results 
from the dimensional reduction. The stationary solu- 
tion can then be written as $ s (0,t) = yj p{9)e l ^ d6v ^ 
and the local dimensionless angular speed of sound 
as c{9) = y/Up{ff)/N . Periodic boundary conditions 
around the ring require the "winding number" w = 
(l/27r) J 2ir dOv(O) to be an integer. 

The qualitative behavior of horizons in this system is 
well represented by the two-parameter family of conden- 
sate densities 

p(0) = ^(1 + b cos 0), 



The mere existence of a black hole solution does not 
necessarily mean that it is physically realizable: it should 
also be stable over sufficiently long time scales. Since 
stability must be checked for perturbations on all wave- 
lengths, the full Bogoliubov [Q spectrum must be deter- 
mined. For large black holes within large, slowly varying 
condensates, this Bogoliubov problem may be solved us- 
ing WKB methods that closely resemble those used for 
solving relativistic field theories in true black hole space- 
times [ Jli| . A detailed adaptation of these methods to 
the Bogoliubov problem will be presented elsewhere p2| . 
The results are qualitatively similar to those we have 
found for black holes in finite traps with low winding 
number, where we have resorted to numerical methods 
because, in these cases, WKB techniques may fail for 
just those modes which threaten to be unstable. 

Our numerical approach for our three-parameter fam- 
ily of black/ white holes in the ring-shaped condensate 
has been to write the Bogoliubov equations in discrete 
Fourier space, and then truncate the resulting infinite- 
dimensional eigenvalue problem. Indeed, writing the 

wave funtion as $ = $> s + ipe 1 -f d8v ^\ decomposing the 
perturbation ip in discrete modes 



<p(9,T) 



uj.n 



(9)+e il 



1 A* v* 



(0), 



and substituting into the Gross-Pitaevskii equation we 
obtain the following equation for the modes and 



where b £ [0, 1]. Continuity, de(pv) — 0, then determines 
the dimensionless flow-velocity field 



v(0) 



UwVl 



27rc((9) 2 



which depends on w as a third discrete independent pa- 
rameter. Requiring that & s (9,t) be a stationary solu- 
tion to Gross-Pitaevskii equation then determines how 
the trapping potential must be modulated as a function 
of 9. All the properties of the condensate, including 
whether and where it has sonic horizons, and whether 
or not they are stable, are thus functions of U, b and w. 
For instance, if we require that the horizons be located at 
9h = ±7r/2, which imposes the relation U = 2ttw 2 (1 — b 2 ), 
then we must have c 2 — v 2 positive for 9 6 (— 7r/2, 7t/2), 
zero at Oh = ±7r/2, and negative otherwise, provided that 
U < 2ttw 2 . The further requirement that perturbations 



to 



E 



h+ f 

ll np Jnp 
~fnp h>np 



In this equation, 



fnp 


2nJ 


fonp 


n 2 



.2jt 



'np 



l 

2Wo 



d Q e -i(n-P)6 



pm -- v > { o)±{c(oy + ~ c{9) 

which, after some lengthy calculations, can be written as 

- U fx — X -A 

h n P = + p) w V 1 ~b 2 a n - p 
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j>|i|, even 

*- £ (^)'(< !+ W ( '' +i) - 

J>\i\, i+J even 

Eliminating Fourier components above a sufficiently high 
cutoff Q has negligible effect on possible instabilities, 
which can be shown to occur at relatively long wave- 
lengths. We face then an eigenvalue problem for the 
2(Q + 1) x 2(Q + 1) matrix built out of blocks of the 
form 

(^np fnp \ 
—fnp h np J 

The numerical solution to this eigenvalue equation, to- 
gether with the normalization condition 
/ d9(ul. <n Uu,', n ' ~ <n v u>' ,n>) = <W<W, provides the 
allowed frequencies. Real negative eigenfrequencies for 
modes of positive norm are always present, which means 
that black hole configurations are energetically unstable, 
as expected. This feature is inherent in supersonic flow, 
since the speed of sound is also the Landau critical veloc- 
ity. In a sufficiently cold and dilute condensate, however, 
the time scale for dissipation may in principle be made 
very long, and so these energetic instabilities need not be 
problematic [ p3[ . 



U 




FIG. 3. Stability diagram for winding number w = 7. Solid 
dark-grey areas represent the regions of stability. Smaller 
plots at higher resolution confirm that the unstable 'fingers' 
are actually smooth and unbroken. Points on the dashed 
curve are states with horizons at Oh = ±7r/2, so that the 
black/white hole fills half the ring. 



More serious are dynamical instabilities, which occur 
for modes with complex eigenfrequencies. Since the Bo- 
goliubov theory is based on a quantized Hamiltonian that 
is Hermitian, there are certainly no complex energy eigen- 
values; but the natural frequencies of normal modes can 
indeed be complex [in which case the usual rule, that en- 
ergy eigenvalues are Ti(n + 1/2) times the mode frequen- 
cies, simply breaks down]. A detailed discussion of the 
quantum mechanics of dynamical instability is presented 
in the Appendix; for the purposes of our main discussion 
it suffices to note that complex (mode) eigenfrequencies 
are indeed genuine physical phenomena, and by no means 
a numerical artifact. For sufficiently high values of the 
cutoff (e.g., Q > 25 in our calculations), the complex 
eigenfrequencies obtained from the truncated eigenvalue 
problem become independent of the cutoff within the nu- 
merical error. The existence and rapidity of dynamical 
instabilities depend sensitively on (M, b, w). For instance, 
see Fig. |^ for a contour plot of the maximum of the abso- 
lute values of the imaginary part of all eigenfrequencies 
for w — 7, showing that the regions of instability are 
long, thin fingers in the (U, b) plane. Figure |J shows the 
size of the largest absolute value of the instabilities for 
each point on the dashed curve of Fig. |^. It illustrates 
the important fact that the size of the imaginary parts, 
which gives the rate of the instabilities, increases start- 
ing from zero, quite rapidly with 6, although they remain 
small as compared with the real parts. 

1 1 1 1 1 1 1 1 1 




FIG. 4. Stability diagram for black/white holes of maxi- 
mum size, i.e., along the dashed line of Fig. ^. 

B. Creation of a black/white hole 

The stability diagram of Fig. |^ suggests a strategy for 
creating a sonic black hole from an initial stable state. 
Within the upper subsonic region, the vertical axis 6 = 
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FIG. 5. Simulation of creation of a stable black/white hole 
and subsequent evolution into an unstable region. Figures 
a-d are snapshots taken at the initial time (a), at an inter- 
mediate time, still within the subsonic region (b), when the 
black/white holes of maximum size is approached (c), and af- 
ter a long time in that configuration (d). Then the parameters 
are changed along the dashed curve of Fig. ^ to enter an un- 
stable region (e) and kept there(f-i). It can be observed that 
a perturbation grows at the black hole horizon and travels 
rightwards until it enters the white hole horizon. 



corresponds to a homogeneous persistent current in a 
ring, which can in principle be created using different 
techniques [p3 . Gradually changing U and b, it is possi- 
ble to move from such an initial state to a black/white 
hole state, along a path lying almost entirely within the 
stable region, and only passing briefly through instabili- 
ties where they are sufficiently small to cause no difficulty. 

Indeed, we have simulated this process of adiabatic cre- 
ation of a sonic black/ white hole by solving numerically 
(using the split operator method) the time-dependent 
Gross-Pitaevskii equation (Q) that provides the evolution 
of the condensate when the parameters of the trapping 
potential change so as to move the condensate state along 
various paths in parameter space. One of these paths is 
shown in Fig. || (light-grey solid line): we start with a 
current at w = 7, b — 0, and sufficiently high U [Fig. 
||(a)]; we then increase b adiabatically keeping U fixed 
until an appropriate value is reached [Fig. ||(b)]; finally, 
keeping b constant, we decrease U adiabatically (which 
can be physically implemented by decreasing the radius 
of the ring trap), until we meet the dashed contour for 
black holes of comfortable size [Fig. 0(c)]. Our simu- 
lations confirm that the small instabilities which briefly 
appear in the process of creation do not disrupt the adi- 
abatic evolution. The final quantum state of the con- 
densate, obtained by this procedure, indeed represents 
a stable black/white hole. We have further checked the 
stability of this final configuration by numerically solv- 
ing the Gross-Pitaevskii equation (Q) for very long pe- 
riods of time (as compared with any characteristic time 
scale of the condensate) and for fixed values of the trap 
parameters. This evolution reflects the fact that no imag- 
inary frequencies are present, as predicted from the mode 
analysis, and that the final state is indeed stationary 
[Fig. p)(d)]. Once the black/white hole has been cre- 
ated, one could further change the parameters (U, b) so 
as to move between the unstable 'fingers' into a stable 
region of higher b (a deeper hole) . 

C. Quasiparticle pair creation 

Instead of navigating the stable region of parameter 
space, one could deliberately enter an unstable region 
[Fig. |(e)-§(i)]. In the this case, the black hole should 
disappear in an explosion of phonons, which may be easy 
to detect experimentally. Such an event might be related 
to the evaporation process suggested for real black holes, 
in the sense that pairs of quasiparticles are created near 
the horizon in both positive and negative energy modes. 
We will explain this point briefly; a more detailed expo- 
sition is included in the Appendix. 

In the language of second quantization, the perturba- 
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tion field operator if satisfies the linear equation 

1 „ . , 



-2 V 



2 c 



1 f 
2 V 



c 2 ] if + c 2 tf>\ 



which, taking into account that \<p{ff), <p' (0')] = 5(8- 
can be written as 

i<p = [<p,H], 
where the Bogoliubov Hamiltonian is 



H = d6 



It// • t / f lc " 

--ipT(f - IViplif + I - — 



i 



-v' + c 2 I (f'tp 



+-(<PV + w ) 



(2) 



The Hermiticity of the Bogoliubov linearized Hamil- 
tonian implies that eigenmodes with complex frequen- 
cies appear always in dual pairs, whose frequencies are 
complex conjugate. In the language of second quantiza- 
tion, the linearized Hamiltonian for each such pair has 
the form 



H 



and the only nonvanishing commutators among these op- 
erators are [A^^, n ,] = 5 nn >. The asterisk on the 
subscript is important: the mode with frequency to* is a 
different mode from the one with frequency w, and A^, n 
is not the Hermitian conjugate of A^ n . It is therefore 
clear that none of these operators is actually a harmonic 
oscillator creation or annihilation operator in the usual 
sense. However, the linear combinations 



A u 







A 



t 



and their Hermitian conjugates are true annhilation and 
creation operators, with the standard commutation rela- 
tions, and in terms of these the Bogoliubov Hamiltonian 
becomes 



H 



E 



[Re(w)(a^a„ - tf n b n ) - Im(a;)(aJ l &£ + a n b n )] 



This interaction obviously leads to self-amplifying cre- 
ation of positive and negative frequency pairs. Evapora- 
tion through an exponentially self-amplifying instability 
is not equivalent, however, to the usual kind of Hawk- 
ing radiation |pH| ; this issue will be discussed in detail 
elsewhere. 



IV. SINK-GENERATED BLACK HOLES 

Condensates which develop black hole behaviors by 
means of flows generated by laser-driven sinks also 



present regions of stability and instability in parameter 
space and, in this sense, their behavior is analogous to 
that in a ring-shaped trapping potential. Here we present 
a simple model that exhibits the main qualitative features 
of more general situations and that can be studied ana- 
lytically. Although in this model, we study a condensate 
of infinite size, in more realistic models or experiments, 
it will suffice to take condensates which are sufficiently 
large, since the stability pattern is not significantly af- 
fected by the (large but finite) size of the condensate. 



A. The model 

Let us consider a tight cigar-shaped condensate of infi- 
nite size such that the motion in the (y, z) plane is effec- 
tively frozen. In appropriate dimcnsionless units, the ef- 
fectively one-dimensional Gross-Pitaevskii equation thus 
becomes: 



id T <P 



1 



d 2 x + v cxt +um 2 $ 



with the normalization condition 



lim 



1 



dx\§(x)\ 2 = n. 



D^oo 2D J_ D 

In this equation, V ex t is the dimensionless effective poten- 
tial that results form the dimensional reduction, which 
already includes the chemical potential. 

In order to obtain a black hole configuration, let us 
choose the potential V cx t so that it produces a profile for 
the speed of sound c(x) = \JlAp[x) of the form 



co, 

c(x) = { cq[1 + (c 

(TCo, 



l)x/e] 



\x\ < L 

L < \x\ < L - 

L + e < \x\ 



with a > 1, and a flow velocity in the inward direction. 
The continuity equation then provides the flow velocity 
profile 



v(x) 



vqc 2 , 



c{xf 



where vq is the absolute value of the flow velocity in the 
inner region. 

As it stand this model fails to fulfill the continuity 
equation at x — 0. In order to take this into account, we 
will also introduce a sink of atoms at x — that takes 
atoms out of the condensate (this can be physically im- 
plemented by means of a laser) . From the mathematical 
point of view, it can be modeled by an additional term 
in the equation of the form —iE5(x) which indeed in- 
duces loss of atoms at x = 0. Equivalently, it can be 
represented by boundary conditions of the form 
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$(0+,t) - $((T,t) = 0, 
$'(0+,t) - $'((T,t) = -2iE$(0,-r), (3) 

which determine the flow velocity inside in terms of the 
characteristics of the outcoupler laser, namely, «o = E. 

Perturbations <fi around this stationary state $ s = 
^p e l J ' v ( x ) dx ^ sucn that $ = $ s + <fi (note that for conve- 
nience we have chosen a different convention as compared 

with the ring in which $ = $ s + ipe 1 J v ) must satisfy the 
boundary conditions (||) and the equation 

i$ = - X -4>" + (c 2 - t> 2 /2 + c"/2c)4> + c 2 e 2l T V* 



where 



a- 1, 



H-Lj-^dil-L-e)]. 



As a further simplifying assumption, we will assume that 
vqc <C 1 so that 



dx'v(x') 



< 



dx 



+ l)x/e] 



< v e< f. 



Let us now expand the perturbation <f> in modes 



Then, the modes u u ^ k {x) and v Ujk {x) satisfy, in each re- 
gion, the Bogoliubov equations 

1 // i/2 2 , 2 2i v 

" / 2 2 /o\ 2 — 2i f X d / j \ 

= 2 U u),fc _ l c -w/2)^,/c-ce J u Wife . (4) 



In the intermediate regions L < \x\ < L + e, the factors 
e ±2i J v ^ ^ c j ag ^ t erms Q £ jrjqg^ jQj become I + 0(e). 
Then the solution of these equations is 

Uoj.k = cvfc + f3uj,kx/e + 0(e 2 ), 

= 7^,fe + n^^x/e + 0{e 2 ) (5) 

as can be easily seen by defining the variable q — x/e so 
that the equations become 



The singular character of c" /c at |x| = L, L + e can 
be substituted by matching conditions at \x\ = L, L + e 
which will in turn provide the connection formulas be- 
tween the modes outside (|x| > L) and the modes inside 
(| a; | < L). Furthermore, the symmetry of the problem 
allows us to study the region x > 0. 

These matching conditions are 

<I>{L + ) — <t>{L~) = 0, 
<£'(L+) - 0'(IT) = — 0(L), 



(f){L + e + )-(i>{L + E-) =0, 
(j>'(L + £+)-</>'(£ + £-) 



cr- 1 



These equations together with the form of the modes in 
the region L < x < L+e, provide the connection formulas 
between the inside and outside modes (from now on we 
will drop the subindex to): 

u- m< k(L) = ~eu' out k (L) + -u out k (L), 

a 

u Lk( L ) = au 'out,k( L ), 

and likewise for the modes v- ln out . 



B. Matching conditions 

The intermediate regions L < \x\ < L + e, provide the 
connection between the perturbation modes in the inner 
and outer regions. Once these connection formulas have 
been established, in the limit of small e, we will only need 
to study the inside and outside modes and their relation 
through such formulas. The case of an abrupt horizon, 
in which the background condensate velocity is steeply 
and linearly ramped within a very short interval, is obvi- 
ously quite special; and it does not particularly resemble 
the horizon of a large black hole in Einsteinian gravity. 
But the connection formula that we derive for this case 
will qualitatively resemble those that are obtained, with 
considerably more technical effort, for smoother horizons 
p2[ . And the results we will obtain for the global Bogoli- 
ubov spectrum of the condensate black hole will indeed 
be representative of more generic cases. 



C. Dispersion relation and boundary conditions for 
large x 

In each of the regions (inside and outside), we can write 

u k (x) = u k e^ k -^ x - L \ v k (x) = v k e^ k+ ^ x ~ L \ 

Upon substitution of this expansion into the Bogoliubov 
equations (Q), we obtain, for each region, the following 
set of algebraic equations: 

h^Uk + c 2 v k = 0, c 2 u k + h^v k = 0, 

where = k 2 /2 + c 2 ± (k\v\ +lu). For these equations to 
have a solution, the determinant must vanish thus pro- 
viding the dispersion relation 

k 4 /4+{c 2 -v 2 )k 2 -2u\v\k-uj 2 = (6) 
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which, for fixed uj, is a fourth-order equation for k. For 
each of the four solutions Uk and Vk must be related by 

Vk = hkUk with hk = 7r(k 2 /2 + c 2 — k\v\ — uj). 

c l 

The constant coefficients can be regarded as normal- 
ization constants and will be set to unity. Let us study 
the possible solutions to the dispersion relation depend- 
ing on whether uj is real or complex. 

Complex frequencies. — In this case all four solutions 
are pure complex, two of them with positive imaginary 
part and two of them with negative one. 

In order to prove this statement, let us first assume 
that there exists a real solution k for a complex u> — 
vo + ij. Then the imaginary part of Eq. (^) implies that 
zu = — k\v\. Introducing this result into its real part 
leads to j 2 — — (k 4 /4 + c 2 k 2 ) which is impossible to fulfill 
because k is real. So the four solutions are complex. 

Because of continuity, all uj in the upper half complex 
uj plane have the same number of solutions with positive 
imaginary part. Otherwise, for some uj there should ex- 
ist a real solution that interpolates between positive and 
negative imaginary part solutions, but this is not possi- 
ble, as we have seen. 

Now let us concentrate on small frequencies, i.e., on fre- 
quencies around uj — 0. For uj — 0, we have a double root 
at k = 0. The other two solutions are k = ±2^fv 2 — c 2 , 
which are real for c 2 < v 2 (i.e., inside) and pure imagi- 
nary for c 2 > v 2 (i.e., outside). Let us follow these four 
solutions when uj = ie. The solutions coming from the 
double root k — will now be of the form k = k r + ieki. 
It is easy to see that k r = and fc^ = l/(— \v\ ± c). If 
c > one is positive and one is negative. If c < \v\, 
both of them are negative. On the other hand, the solu- 
tions k = ±2Vw 2 — c 2 , for c 2 > v 2 , are already complex 
conjugate. For c 2 < v 2 , we write k = k r + ieki and in- 
troduce it into Eq. (||). We then see that at first order 
in e, ki = \v\/(v 2 — c 2 ) > 0. Thus, we have seen that for 
uj = ie, we have two solutions with positive imaginary 
part and two with negative one in any case (inside and 
outside). But if this is so for uj = ie it must be true in 
the whole upper uj plane and consequently in the whole 
complex uj plane. □ 

In the inside region all possible solutions are in prin- 
ciple allowed but outside we are only left with the two 
that have Im(fc) > 0, because the other two grow expo- 
nentially. 

Real frequencies. — Outside (c 2 > v 2 ), there are two 
real and two complex conjugate solutions. Of these two 
complex solutions, only one is allowed [the one with 
Im(fc) > 0] because the other grows exponentially. In- 
side (c 2 < v 2 ), for uj > w max there are two real and two 
complex conjugate solutions; for uj < w max there are four 
real solutions; the value uj = cj max is a bifurcating point. 



D. Connection formulas for complex frequencies 

Since we are interested in the existence of dynamical 
instabilities, we will concentrate in the case in which uj 
is complex. Then, as we have seen, the dispersion equa- 
tion ([|) has four complex solutions for k in each region. 
Inside, all four solutions fci n .i, i = 1---4 are in prin- 
ciple possible but outside those with Im(fe out ) < will 
increase exponentially. Therefore, up to corrections com- 
ing from the finite size of the condensate, which we ig- 
nore here, only modes associated with /c ut,a, ct = 1,2 
such that Im(fcout,a) > are allowed. Each mode 
u ou t,a(x) — e l< - kout ' a ~ v °/' j2 " >i ' x ~ L " > will match a linear com- 
bination u- m , a (x) — J2i FaiUin,i{x) of modes u- lai i(x) in- 
side, i.e., 

ii- (r) — F -P l ( k ^.*- V a)( x - L ) 
"an, a y**} — / ± at*-- 

i 

and similarly for w ou t,Q and v- mta : 

v (r) = V* F h- A k ^,i+ v o)( x - L ) 

i 

After some straightforward calculations, it can be seen 
that these connecting coefficients F a i are given by F a i = 
Y,j(M- 1 ) ij C aJi where 



1 


1 


1 


1 


in,l 




^in,3 


^in,4 


in, 1 


hin,2 


^in,3 


hin,4 



I l/a - iek~ nt<a \ 

\ ^out.a^out.a / 

In these equations, 

fc in.i = foiM ± «o and fc^ ut:Q = fc out:Q ± vq/o- 2 . 

E. Boundary conditions at the sink. Complex 
eigenvalues 

We have already found the modes in the inner and 
outer region as well as their relation. To determine which 
(complex- frequency) modes will be present, there only 
remains to impose the boundary conditions dictated by 
the presence of the sink at x = 0. 

As we have already mentioned, the symmetry of the 
system under reflection (x — > —x) allows us to study 
only the region x > 0, provided that we study the even 
and odd perturbations separately. For odd fluctuations 
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come 



, t)], the boundary conditions (^) be- where fco = — c\ . More explicitly, for L smaller than 



0o(O,r) =0 



at all times r. This implies that the u and v compo- 
nents of 4> must separately satisfy the boundary condi- 
tion. Since we can have any linear combination of the 
two solutions that decay outside the horizon, we there- 
fore have a two-by-two matrix constraint. The condition 
that a non-zero solution exists is that the determinant 



det 



"in, l(0) Uin : 2(0) 
fin, l(0) Win, 2(0) 







and therefore 

^2F li F 2j (h ia , i - h iaJ )e- i ^+ k ^ L = . (7) 
ij 

For even fluctuations [4> c (x,t) — e (— x, r)], the 
boundary conditions (|^) become 

&(O,t)+™ o 0c(O,t) = 0, 



which implies that 



-i(k iTL ,i+ki n j)L 



0. (8) 



For fixed L, U, Vq, and cr, the quantities F, /ii n , and 
fein that appear in Eqs. ([?]) and (|) are only functions 
of u>. Therefore the solutions to these equations are all 
the possible complex eigenfrequencies, which depend on 
the free parameters that determine the model, namely, 
the size 2L of the inner region, the speed of sound inside 
Co, the relative change of the speed of sound between 
the inner and the outer regions cr, and the flow velocity 
inside Vq (related to the characteristics of the outcouplcr 
laser). In practice, there are also other parameters of 
the condensate such as its size 2D (which has been made 
arbitrarily large) and the size of the intermediate regions 
e (which has been made arbitrarily small) . 

Equations @) and (||) can be solved numerically for 
different values of the parameters cr, U, vq, and L. The 
numerical method employed is the following. The equa- 
tions above have the form 

f(w;a,U,v ,L) = , 

where / and oj are both in general complex. We plot 
contours of constant absolute value of / in the complex 
lu plane; where |/| approaches zero, we have an eigenfre- 
quency. 

The distribution of complex solutions in the complex 
lu plane depends on the size of the inner region L, for 
given cr, U, and vq. Direct inspection of the numerical 
results shows that the number of instabilities increases 
by one when the black-hole size L is increased by 7r/fco 



7t/ko — 5 (5 being much smaller than Tr/ko) there are no 
complex eigenfrequencies; for (L + <5)fco/7r € [n,n + 1] 
with n = 1,2..., we have n complex solutions except 
for L = (n + l/2)7r/fc , where we find n — 1 complex 
solutions instead of n [i.e, there is one mode for which 
Im(w) = within numerical resolution]. This can be eas- 
ily interpreted qualitatively since the unstable modes are 
basically the bound states in the black hole, and the high- 
est wave number k on the positive norm upper branch, 
for the barely bound state with lu — > _ , is exactly fco. 
So the threshold is simply when the well becomes big 
enough to have a bound state; the small 5 displacement 
comes in because the horizon is not exactly a hard wall; 
and similarly for the extra bound state every 7r/fco. Thus 
stability can only be achieved for small sizes of the in- 
ner region, L < 7r/fco. As we discussed in Sec. II, the 
wave length 27r/fc of the perturbations must be smaller 
than this size, which implies fc > 2ir/L > 2fco. How- 
ever, for these perturbations the hydrodynamic approx- 
imation, which requires fc < 2fco, is not valid. Therefore 
there are no stable black hole configurations in a strict 
sense. The sizes of the imaginary parts of the complex 
solutions decrease as the size L of the interior of the black 
hole increases. Thus although a larger hole has more un- 
stable modes, it is actually less unstable (and might even 
became quasistable in the sense that its instability time 
scale would be longer than the experimental duration). 



V. CONCLUSIONS 

We have seen that dilute Bose-Einstein condensates 
admit, under appropriate conditions, configurations that 
closely ressemble gravitational black holes. We have an- 
alyzed in detail the case of a condensate in a ring trap, 
and proposed a realistic scheme for adiabatically creat- 
ing stable sonic black/white holes and we have seen that 
there exist stable and unstable black hole configurations. 
We have also studied a model for a sink-generated sonic 
black hole in an infinite one-dimensional condensate. The 
dynamical instabilities can be interpreted as coming from 
quasiparticle pair creation, as in the well-known sug- 
gested mechanism for black hole evaporation. General- 
izations to spherical or quasi-two-dimensional traps, with 
flows generated by laser-driven atom sinks, should also be 
possible, and should behave similarly. While our analysis 
has been limited to Bogoliubov theory, the further the- 
oretical problems of back reaction and other corrections 
to simple mean field theory should be more tractable for 
condensates than for other systems analogous to black 
holes. And we expect that experiments along the lines we 
have proposed, including both creation and evaporation 
of sonic black holes, can be performed with state-of-the- 
art or planned technology. 
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APPENDIX: COMPLEX FREQUENCIES. 
REDUNDANCY AND NORMALIZATION 

In this Appendix, we will analyze the issue of the re- 
dundancy and normalization of the Bogoliubov modes in 
the presence of complex frequencies from a general point 
of view. Dynamical instabilities in quantum field theory, 
and the quantization of dynamically unstable modes, do 
not seem to be widely understood: it is for instance com- 
mon to read axiomatic statements that one must only 
quantize positive norm modes, even though this implic- 
itly neglects dynamical instabilities, and does not follow 
in general from the fundamental commutation relations. 
But some explicit treatments of quantum instabilities 
have been available in the literature for some time [p5[ ; 
here we review this subject in the specific Bogoliubov 
context. 

We will begin by writing the Bogoliubov equations in 
their most usual form 



TlLU.; 



/i (x) c(x) 2 e 2W « 
-c(x) 2 e 2l,5 ( x ) -ho(x) 



(Al) 



where h (x.) — —h 2 /2mV 2 + U cxt (x) + 2mc(x) 2 — fi. In 
terms of these modes, the atomic second quantized field 
operator has the well-known form ^(x, t) = ^(x, t) + 
■0(x, t) with 



<Kx,i) = [a^(x) e -^* + aj« 3 -(x)* 



At 



(A2) 



If there is a solution (uj,Vj) to Eq. ( Al ) with mode fre- 



quency ojj , then straightforward substitution shows that 
(Uf,Vj>) — (Vj,Uj) must be a solution with frequency 
u>ji = —lj* . If we examine the contributions of these two 
solutions, however, we find that together they yield but 
a single term, of the form [dj + at,)e~""*itj(x) + (at + 
j*t„ 



Ij, 



)e v 



(x). It is thus a quite trivial fact that the two 



modes j and j' are redundant. We are free to simplify 



our notation by redefining dj 



an, and eliminat- 



ing mode j' (leaving it out of the sum over frequencies) . 
Alternatively we could of course eliminate j and keep j'. 
Which of these two notational conventions we should take 
is best determined by the commutator [dj + at, , a] + dj>], 
which will tell us whether the coefficient of Uj — v*i is 
properly an annihilation operator or a creation operator. 



Since the only commutation relations that we are given 
are those of ip and we must derive the orthogonality 
relation for solutions of Eq. (Al), and use it to invert 



Eq. (A2). We can use Eq. (Al) to show that 



(u>j + u>k)Mjk = I d 3 x (ujV k ~ VjU k ) = 0, 

3, 



{u 3 - uj* k )N jk = (luj - w%) J d x (ujU* k - Vjv k ) = , (A3) 

where in the case of infinite volume the right hand 
sides are zero in the distributional sense, being infinitely 
rapidly oscillating boundary terms. This obviously im- 



plies that Mj k vanishes unless iv k 



and Njk van- 



ishes unless uj k = u>j . One can then show that it is always 
possible to take linear combinations among degenerate 
modes, and to eliminate redundant modes as just dis- 
cussed, in such a way as to make Mj k always vanish, and 
Njk = S k j, where for every j there is a single dual mode 
j, with luj = lu* . In the case of real Uj , but only then, we 
have j = j. In general, though, duality is reciprocal (the 
dual mode of j is always j). 

The result is that we can now insert Eq. (A2) into 
the second quantized Hamiltonian, with the T-matrix 
approximation for the interparticle interaction, to obtain 
the linearized Bogoliubov Hamiltonian for the perturba- 
tions 



(A4) 



Since the sum over all modes j also includes the dual to 
every mode with complex ojj , H is manifestly Hermitian 
even though u>j need not be real. We can also invert Eq. 
3) to learn that 



d x [ujS'ijj + VjS'ip 



(A5) 



which with Eq. (A3) implies the commutation relations 



[aj,a k \ 







(A6) 



For all j with real Wj , Eq. ( A6 ) are merely the standard 
canonical commutation relations; and our normalization 
conventions Mj k — and Njk — Sj k are likewise the ones 
most often presente d. I n the case of complex tOj where 
j^j, however, Eq. ( |A6| ) implies that the canonical con- 
jugate of dj is al, and this is no longer the same as the 
Hermitian conjugate at. In fact for complex u>j we have 
= 0; this already follows from the second line of 



Eq. (A3), which implies that the norm Njj of any mode 



with complex ujj is zero. But if dj and aj commute, then 
it is clear that neither dj nor dj is really a harmonic os- 
cillator annihilation operator in the usual sense, nor are 
d\ or d\ proper creation operators. The commutation 



relations (A6) are validly derived from the fundamental 
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commutation relations for if; and ifi^; but they do not im- 
ply, for instance, that either atfij or djolj has the discrete, 
equally spaced spectrum, bounded from below, that one 
expects of a quasiparticle number operator. 

To understand the dual pairs of modes with complex 
frequencies, we can define the ordinary annihilation op- 
erators 



77= (Sj +%), 



bj- 



1 M 
\/2 V 3 



4) 



(A7) 



and their Hermitian conjugates, among which the only 
nonvanishing commutators are the ordinary 



[b j ,b\] = [b J M = i 



(A8) 



In terms of these operators, which are harmonic oscillator 
annihilation and creation operators with all the familiar 
properties of such, the j,j subsector of the Bogoliubov 
Hamiltonian H appears as 



H j ,j=Re(u; j )[btb j 



Imiujj^Pbl + bjbj] . (A9) 



Note that Eq. (A9) is only the simplest form in which 
one may write the j, J sector of the Hamiltonian: by in- 
troducing appropriate factors of e ±M / 2 / cos a into Eq. 
(A7), for any a, we can make Im(wj ) — > Im(wj)/ cos a 
and add a term Im(u>j) t&na(b^bj + bjbj). 

We can now examine the spectrum of Hjj by consid- 
ering it in the basis of eigenstates of h = b^bj + bib 



and A = Mbj 
defining 



\Ea) 



OO 



3 U 3 



In fact A commutes with Hi so 



A)|r 



A > 



\E A ) = J2 c n\n)\n + A) A < 



(A10) 



n=0 



where b^b^m^n) = m\m)\n) and bjbj\m)\n) = n\m)\n), 



we find Hjj\E A ) = [^ARe(^) 
the recursion relation 



Ea]\Ea}- And we have 



EaCu = lm(uij)[^/n(n + A)cn 



+ + !)(« + A + l)c„ H 



(All) 



As n — > oo, we have c n +i — > — c„_i, and so J2 n |c„| 2 does 
not converge: none of the eigenstates of Hjj is normal- 
izable. One can however obtain delta- function normal- 
ization for a continuous spectrum of real E A , bounded 
neither above nor below. 

That the Hamiltonian H is unbounded from below 
does not indicate anything unphysical about our model: 
we have simply linearized about an unstable excited state 
of the nonlinear full Hamiltonian, which is bounded from 
below. Real negative frequencies ojj, where our conven- 
tion Njj = 1 has been imposed, indicate energetic insta- 
bilities, whereby the system will decay in the presence of 



dissipation. Complex loj, on the other hand, indicate dy- 
namical instabilities. Classically, a dynamically unstable 
system will exponentially diverge from the initial station- 
ary state if is perturbed, even without dissipation. Quan- 
tum mechanically, we have just seen that a dynamically 
unstable system has no normalizable stationary states. 
If an initially stable system is driven into a state which 
is stationary but dynamically unstable at the classical 
(mean field) level, the initial state will have had finite 
Hilbert space norm, and hence under unitary evolution 
the final state will have the same norm. Thus it will not 
be a stationary state; one may say that quantum fluctu- 
ations will always trigger the dynamical instability For 
a logarithmically long period of time, however, the lin- 
earized theory will still remain valid. In this sense, our 
linearized description of quantum dynamical instabilities 
is sound. 
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